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Limits on WWZ and WWj couplings 
from WW and WZ production 
in pp collisions at yfs = 1.8 TeV 

Abstract 

Direct limits are set on WWZ and WWj three-boson couplings in a search 
for WW and WZ production with high transverse momentum in pp colli- 
sions at y/s = 1.8 TeV, using the Collider Detector at Fermilab. The results 
are in agreement with the SU(2) X U(l) model of electroweak interactions. 
Assuming Standard Model WWj couplings, the limits are interpreted as di- 
rect evidence for a non-zero WWZ coupling at subprocess energies near 500 
GeV. Alternatively, assuming identical WWZ and WWj couplings, bounds 
-0.11 < k < 2.27 and -0.81 < A < 0.84 are obtained at 95% C.L. for a form 
factor scale 1000 GeV. 
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Since formulation of the SU(2)xU(l) gauge theory of electroweak interactions [1], 
many of its predictions have been confirmed, including the existence of the W and 
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Z force carriers. The parameters of the model have been determined with ever in- 
creasing precision, but only now are there direct tests of the predicted interactions 
of W, Z } and 7 bosons with each other. These interactions are the most charac- 
teristic and fundamental signatures of non-Abelian symmetry in the theory. The 
predicted interactions are described by trilinear couplings WW Z and WW~/ } which 
we address here, as well as quadrilinear couplings. Boson pair production is sensi- 
tive to these couplings [2], and the WW^f coupling has been tested directly in the 
process pp — > W~/[3]. In this letter we report direct information about the WW Z 
and WW^f couplings obtained from a search for WW and WZ production with large 
boson transverse momentum (Pt) in PP collisions at yfs = 1.8 TeV. 

Indirect limits on the WW Z and WW^f couplings have previously been set based 
on one-loop effects at low energies and precision measurements at the Z resonance[4, 
5]. The direct measurements provided by diboson production are valuable because 
they are unambiguous as to their interpretation. They do not require additional 
theoretical assumptions or calculation of loop diagrams, which can present theoretical 
ambiguities for non-standard models. Furthermore, they are sensitive to directions 
in the space of couplings which are not well constrained by the indirect limits. 

The most general WW Z and WW^f couplings consistent with Lorentz invari- 
ance have been formulated and may be parametrized in terms of fourteen indepen- 
dent couplings (or form factors), seven for the WW Z vertex and seven for the WW^f 
vertex [6]. They are usually denoted g\ , g\ ^ , \ v , k v , \ v , and k v where V is ei- 
ther Z (for WW Z) or 7 (for WW^f). The standard SU(2) X U(l) electroweak theory 
corresponds to the choice gj = g^ = 1 and k 1 = k z = 1 with all other couplings set 
to zero. 

In the Standard Model, the dominant contribution to diboson (WW or WZ) 
production in pp collisions at y/s = 1.8 TeV comes from two types of Feynman di- 
agrams, the t- or u-channel diagrams, which involve the couplings of W and Z to 
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fermions, and the s-channel diagrams which are the only ones containing the three- 
boson coupling. To the extent that the fermionic couplings of the W and Z have been 
well tested, we may regard diboson production as primarily a test of the three-boson 
couplings. There are substantial cancellations between the s-channel and the t- or 
u-channel diagrams, resulting in cross sections of 9.5 pb and 2.5 pb for WW and WZ 
production respectively. If any of the three-boson couplings differ from the standard 
model values then the cancellations are reduced and the cross section increases. The 
enhancement is greatest at high boson Py where the strongest cancellations occur 
in the standard model. Therefore, this analysis looks for anomalously large cross 
sections at high boson Py in order to obtain information on the couplings. 

The data for the analysis were recorded with the Collider Detector at Fermilab 
during the 1992-93 Fermilab Tevatron collider run, corresponding to an integrated 
luminosity of 19.6 pb -1 . The detector has been described in detail elsewhere[7]. Here 
we give a brief description of the components relevant to this analysis. The location of 
the event vertex is measured along the beam direction with a time projection cham- 
ber (VTX). The momenta of charged particles are measured in the central tracking 
chamber (CTC), which is surrounded by a 1.4 T superconducting solenoidal mag- 
net. Outside the CTC, the calorimeter is organized in electromagnetic (EM) and 
hadronic (HAD) compartments with projective towers covering the pseudorapidity 
range \r]\ < 3.6. Outside the central calorimeter, the region \r]\ < 1.0 is instrumented 
with drift chambers for muon identification. 

Each electron is identified by an isolated cluster in either the central EM 
calorimeter (\r]\ < 1.1) which matches a track in the CTC or the endplug EM calorime- 
ter (1.1 < \r]\ < 2.4) with associated hits in the VTX. Each muon is identified by an 
isolated track in the CTC with minimum ionizing energy in the calorimeter. Events 
with one or more muons must have at least one muon with matching hits in the muon 
chambers. The presence of neutrinos is inferred from missing transverse energy ($ T ), 
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which is measured by the magnitude of the vector sum of the calorimeter tower en- 
ergies perpendicular to the beam axis. Jet energy is measured by clustering the EM 
and HAD calorimeter energy within a cone AR < 0.4, where AR = ^A(f) 2 + Ar] 2 , 
and (f) is the azimuthal angle [8]. 

We search for WW and W Z event candidates consistent with the decay of 
one boson to leptons and the other to hadrons. This choice of decay channels gives 
better sensitivity to anomalous three-boson couplings than the purely leptonic chan- 
nels because the leptonic branching fractions of the W and Z are small and because 
the acceptance of the detector for jets is larger than for leptons. Background from 
the QCD processes pp — > W + jets and pp — > Z + jets is greatly reduced by re- 
quiring a large boson Py, while retaining good sensitivity to anomalous three-boson 
couplings[6]. Background QCD processes are calculated at Born level [9], including 
simulation of the CDF detector and jet fragmentation using an adaptation of the 
HERWIG program[10, 11]. The boson Py requirement for WW and WZ event se- 
lection is chosen so that less than one background event is expected in the final 
sample. With this choice it is unnecessary to perform a background subtraction and 
any theoretical uncertainty in the background calculation is avoided. 

A leptonic W decay is identified by an isolated electron or muon with Py > 
20GeV/c and ^ T > 20 GeV forming a transverse mass My > 40GeV/c 2 . A leptonic 
Z decay is identified by an electron or muon pair of opposite charge forming an 
invariant mass 70 < M < 110GeV/c 2 . In events with a leptonic W or Z decay, a 
candidate hadronic W or Z decay is defined by the two jets (leading jets) in the event 
with the highest jet transverse energies (Et)- Each jet must have E? > 30 GeV and 
the invariant mass of the jet pair must be in the range 60 < Mjj < 110GeV/c 2 . 
The Py of the two-jet system, interpreted as a hadronic W or Z decay, is required 
to satisfy Py > 130GeV/c for leptonic W events or Py > lOOGeV/c for leptonic Z 
events. 
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The two-jet mass spectrum is shown in Fig. la for events with a leptonic W 
decay and with both leading jets satisfying E? > 30 GeV. The sum of the predicted 
Standard Model WW and WZ signals plus QCD background is also shown, where the 
background is normalized to the observed number of W events with two jets minus 
the predicted signal. Fig. lb shows the two-jet Py distribution in the subset of events 
which satisfy the two-jet mass criterion. The two-jet Py requirement is indicated by 
the arrow. One event passes this cut. For events with a leptonic Z decay there are 
no events which satisfy all selection criteria. 

The limits on the couplings follow from a Monte Carlo calculation of expected 
event yields for various values of the couplings. The Monte Carlo event generator[6, 
12] calculates to leading order the processes pp — > W + W~ and pp — > WZ with 
subsequent decay of a W to eu, fj>v } or jj and a Z to ee, /i/i, or jj. Higher order QCD 
corrections to the cross section are accounted for by a "K-factor" of K = 1 + |7ra s [6]. 
MTB2 structure functions are used[13]. Initial and final state QCD radiation effects 
and jet fragmentation are modelled with an adaptation of HERWIG[10, 11]. The event 
generator is combined with a detector simulation which includes trigger efficiencies, 
lepton identification efficiencies, and jet response modeling. A fast parametrization of 
the full detector simulation was also employed. The trigger and lepton identification 
efficiencies are determined from the data and amount to 78% for electrons and 79% 
for muons. The modeling of the jet response and resolution are tuned to agree with 
studies of collider and test beam data[14]. The two-jet mass resolution is expected 
to be 9 GeV/c 2 for diboson events that would pass our candidate selection criteria. 
The efficiency of the two-jet mass cut is 88% for events passing all other cuts. 

The systematic uncertainties on the yield are the uncertainties in the structure 
functions (6%), jet E? scale and resolution (16%), luminosity (4%), lepton identi- 
fication efficiency (1%), and trigger efficiency (1%). The Monte Carlo acceptance 
modeling has 3% statistical uncertainty, and a 5% systematic uncertainty allows for 
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differences between fast and full detector simulations. In addition a 14% uncertainty 
is assigned for the effects of higher order QCD corrections[10, 15, 16]. These uncer- 
tainties are combined in quadrature. 

The acceptance is a strong function of the couplings, because of the boson 
Pt cut in combination with a varying boson Py distribution. For standard model 
couplings, 0.13 WW/WZ — > lujj events and 0.02 WZ — > 11 jj events are expected to 
pass the selection criteria, where I is either an electron or a muon. The observation 
of one event in the lujj channel and zero events in the 11 jj channel is therefore not 
indicative of a departure from standard model couplings, even without consideration 
of the QCD background. 

The predicted yield of high Py boson pairs is a quadratic function of the anoma- 
lous couplings. The lack of an excess of events therefore results in bounds on the cou- 
plings which take the form of ellipses in the plane of any two couplings. Since the one 
event passing all selection criteria could be either signal or background, we calculate 
the confidence limits from the probability of observing one or less signal events. We 
do not perform a background subtraction and therefore obtain conservative limits. 
The probability distribution used is the convolution of a Poisson distribution with a 
Gaussian, where the Gaussian smears the mean of the Poisson distribution around 
the expected yield within the systematic uncertainty. 

In the calculation of WW and W Z cross sections, the anomalous parts of the 
couplings are suppressed at high subprocess center of mass energy by a dipole 

form factor [6]: 



where £ stands for any of the couplings, ^sm is its value in the standard model, and 
Aff represents the energy scale of unknown phenomena. Without this suppression, 




e(o) - 
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the anomalous couplings would result in cross sections that violate unitarity at large s. 
With the suppression, the couplings approach their standard model values at energies 
above the scale Aff, and the cross sections respect unitarity as long as the anomalous 
couplings are not too large [17]. 

In Fig. 2 we present bounds on four pairs of couplings. Except as noted in 
the figure caption, for each case all the other couplings are fixed at the standard 
model values. Each pair is constrained to the interior of an ellipse, which is a two 
dimensional section through an ellipsoidal allowed region in the fourteen dimensional 
space of three boson couplings. Because the bosons are required to have high Pj 
our search is most sensitive to the couplings at energies near a/I = 500 GeV. The 
limit contours, however, correspond to the value of the couplings at a/I = and 
therefore depend on the choice of Aff according to equation (1). The bounds are 
shown for Aff = 1000 GeV and Aff = 1500 GeV. The unitarity bounds, which 
depend strongly on Aff, are also shown[17, 18]. For values of Aff larger than about 
1600 GeV the bounds from unitarity are stronger than the bounds from the search. 

Fig. 2a shows limits in the plane A 7 vs. \ z . The limits are stronger for \ z , 
illustrating the fact that the search is in general more sensitive to the WW Z couplings. 
It is therefore complementary to studies of the process pp — > W~/[3]. 

The limits of Fig. 2b focus on the WW Z vertex, assuming that the WW^f 
couplings take their standard model values. Bounds are shown for the couplings g z 
and k z , which are the only WW Z couplings predicted to be nonzero in the standard 
model. The fact that the point g z = k z = lies outside the allowed region can be 
interpreted as direct evidence for a non-zero WWZ coupling, and for the resulting 
destructive interference between s-channel and t- or u-channel diagrams which takes 
place in the standard model. Specifically, the search is directly sensitive to the WWZ 
coupling in the region a/I = 500 GeV. If the WWZ coupling were zero in this region, 
the s-channel diagram containing the WWZ vertex would not contribute to the am- 
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plitude, and the other diagrams by themselves would predict the observation of 15 ±3 
events, where the uncertainty is systematic. Independent of the choice of A FF} this 
possibility is excluded at greater than 99% CL. 

Figs. 2c and 2d show limits on the couplings k and A, assuming specific relations 
between the WW Z and WW^f couplings. In Fig. 2c, the WW Z couplings are assumed 
to equal the WW^f couplings. The resulting 95% CL limits on k and A separately, 
assuming that only one departs from its standard model value, are —0.11 < n < 2.27 
and -0.81 < A < 0.84 for the choice A FF = 1000 GeV. With the assumption of 
matching WW Z and WW^f couplings, limits also result for the W boson electric 
quadrupole moment QY = M " e 2 (k — A) and magnetic dipole moment /i w = 2 ^ (1 + 
k + A). In the standard model, these moments take the values QY = M " e 2 and 
ld W = jj^- The point QY = = is outside the allowed region. Assuming only 
one of the moments departs from its standard model value, the limits at 95% CL are 
-2.42 < QY/{e/M w 2 ) < 0.35 and 0.37 < /j w /(e/M w ) < 1.70 for A FF = 1000 GeV. 

For Fig. 2d, the relation assumed between the WW Z and WW^f couplings is 
given by the HISZ equations[4], which specify \ z , k z , and g z in terms of the in- 
dependent variables n 1 and A 7 . This prescription preserves SU(2) X U(l) gauge 
invariance and is well motivated in an effective Lagrangian approach. The corre- 
sponding subspace of anomalous couplings is not well constrained by previous indirect 
measurements [4]. The individual 95% CL bounds on A 7 and n 1 are —0.35 < n 1 < 2.57 
and —0.85 < A 7 < 0.81 for A^ = 1000 GeV, if only one of the two is varied from its 
standard model value. 
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Figure 1: Selection of WW/WZ — > lujj candidates. All event selection cuts except 
the two-jet mass and two-jet Py cuts were used to select the events in (a). The 
subset of events from (a) passing the two-jet mass cut is shown in (b). One event 
remains after making all cuts. The solid line shows the data, the dots show the 
predicted standard model diboson signal, and the dashes show the predicted signal 
plus background shape. 
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Figure 2: Allowed regions for pairs of anomalous couplings. All couplings, other than 
those listed for each contour, are held at their standard model values. The solid lines 
are the 95% CL limits and the dotted lines are the unitarity limits; each is shown for 
Aff = 1000 GeV (outer) and 1500 GeV (inner). The + signs indicate the Standard 
Model values of the couplings, (a) A 7 and \ z ; (b) g\ Z and k z ; (c) k and A assuming 
the WW Z and WW^f couplings are the same; (d) k 7 , k z , A 7 , \ z and g z in the HISZ 
prescription (see text), with independent variables n 1 and A 7 . 
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